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The use of the classical path approximation in line broadening theory is illustrated in a develop- 
ment of the familiar impact and one-electron theories. The one-electron theory which is presented 
is an improved version of the usual one-electron theory: this improved version of the theory provides 
a consistent description of a line profile from the half width to the quasi-static wings. The validity criteria 
for the impact and one-electron theories are discussed in detail and a comparison of the theories is 
made with a view toward the development of a more general theory containing elements of both. 
To avoid unnecessary mathematical complications and to provide a more transparent comparison 
of the theories, the Stark broadening of the Lyman series in hydrogen is used as a specific example. 
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1. Introduction 

In a previous paper [l], 1 henceforth referred to as 
I, a general expression for the thermal average, eq 
(1—47), was derived using elassieal path methods. In 
this paper we will show how this expression may be 
used to derive the familiar results of the impact 
theories [2, 3, 4] and the one-eleetron theories [3, 5]. 
Possible improvements to these theories by means of 
a more general treatment of the thermal average, will 
also be discussed. 

The perturber coordinates, x (eigenvectors of the 
perturber position operator) are integrated out of the 
expression for the thermal average when classical 
path methods are used (see eq (1-36)). The only per- 
turber variables which need to be considered in the 
expression for the thermal average, eq (1-47) are the 
coordinates x (/) , v (/) of the perturber wave packets 
(which are viewed as classical particles). The super- 
script {i) notation, which was used in I to distinguish 
between these variables, is no longer needed and. for 
convenience in notation, this superscript will not be 
used in this paper. 

The thermal average will be evaluated by means ol 
statistical techniques. In order to clearly illustrate 
these techniques and avoid unnecessary mathematical 
complications, we will consider the case of the Lyman 
lines emitted by hydrogen atoms in a plasma. 
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2. Preliminary Mathematics and 
Approximations 

2.1. The Line Shape in the Static Ion Approximation 

We will make the usual quasi-static approximation 

[2, 3] for the ions in the plasma. The electric field, ST, 
produced by these ions is regarded as being static 
during the time of interest, and this static field is used 
to define the Z axis for the atom. H n is the Hamiltonian 
for an isolated hydrogen atom and, for convenience, 
the perturbation is taken to be a dipole interaction 
with the atomic dipole — eR, 



V(t)=eZV+V e (Ri x, v, /) 
F,.(R, x, v, t) = eR -gV(x, v, t) 



(1) 

(2) 



where Z denotes the Z component of R, and % e (t) is 
the fluctuating electric field produced by the elec- 
trons. The classical variables x and v are 3/V-vectors 
describing the N electron perturbers (these were 
denoted by x (,) and v (/) in I). The average over ion 
holds is performed by means of the usual microfield 
average [2, 3, 6]. In this case, the complete line profile 
is given by 



■J 



Xo>)= ^(gr)/(oi,gr)r/gr, 



where S^^S) is the probability of finding 
of magnitude 8f. The function 7(a>, jf) is 



(3) 

an ion field 
the same as 
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the function I{oj) discussed in I, except that we now 
explicitly indicate a dependence on the ion field £?. 
That is, we may regard the radiator as being an atom 
which is subjected to a static field £?: the perturbers 
are electrons and the resulting line shape, for a given 
fields, is /(co,£f): the total line profile is then obtained 
by averaging /(a>, gO over all possible ion fields as 
shown in eq (3). 

2.2. The Time Development Operator in the No- 
Quenching Approximation 

The time development operator T ( ,(t) (see eq (1-39)) 
is the solution of the differential equation 



ih — T„(t) = [H n + eZZ + VAt)]T n {t) (4) 

at 



T„(0) = l. 
In an interaction representation defined by 

T„(t) = exp(-itH„lh)T;,(t), (5) 

we have 

ih^T;,(t) = [egz'(t) + v;u)]T;,(t) (6) 

at 
Z'U) = exp {itH (l lh)Zexp(-itHalfi) (7) 



state neglecting spin) which we shall denote by |1 > ; 
the frequency of the line of interest may then be de- 
noted by co,, i = (E„ — E\ )jh. 

A radiationless transition which depopulates the 
initial level E „ may be said to quench the radiation 
near co,, l.The probability for such an inelastic transi- 
tion from an initial state \///> to some other state 
\/t'Ji' > (where /* t^ // ) is given by the square of 
(/t/'\T n | //A'). The no-quenching approximation states 
that these "off-diagonal" (i.e., n ^ /*' ) terms may be 
neglected (see eq (12) of ref. |9|). The justification for 
this approximation will be given in section 7; in this 
section we will simply use <///-\T,,\///-' > = to 
simplify the form of T„(t). 

Since the ' 'off-diagonal" matrix elements of T„{t) 
are assumed to be negligible, it is clear from eq (5), 
that the "off-diagonal" elements of 77,(0 must also 
be negligible. Noting that the H a eigenvalues depend 
only on /• . the differential equation for 77,(0 . eq (6), 
becomes 

ih{-{,J\r„(t)\,J') 

at 

= ^i"A<'%z + v,,(t)\,j")(, l / , '\r ll (t)\,,/.') (io) 

with the aid of eqs (7) and (8). If we define a projection 
operator P which picks out the part of any operator 
which is diagonal in n 

(,J\PZ\„'/')=8„„'(,J\Z\,</'), 

then eq (10) may be written 



17,(0 = exp (itH„lh)V e (t) exp(-itH„lh). 
The formal solution of eq (6) may be given by 

Tl(t) = 0exp\-(M) i[e$Z'(t') + VUt')]df 



(8) 



(9) 



where 6 is the time ordering operator [7] which is 
required since [egfZ'(0 + 1^(03 does not commute 
with \_e%Z'{t') + V' e {t')] unless t = t' (cf. eq (6) of 
ref. [8]). Equation (9) is a very complicated expression 
and approximations must be employed to simplify the 
form of this operator. In this section we will consider 
the no-quenching approximation [9] which is frequently 
used for hydrogen lines. 

To make an explicit statement of the no-quenching 
approximation, we must first introduce a complete 
set of//,, eigenstates. The H„ eigenvalues,/?,,, depend 
only on the principal quantum number//, hence we 
denote the f/„ eigenstates by | ••/> where y£ is some set 
of quantum numbers which index the degenerate states 
of the level E„. We are interested in the radiation pro- 
duced when the atom makes a radiative transition 
from a group of initial states (degenerate states of 
some particular level E„) to a group of final states 
which have a lower energy. For the Lyman series in 
hydrogen, there is only one final state (i.e., the ground 



ih-T„{t) 

at 



[e%PZ + PV e {t)]T' n {t). 



(11) 



The appearance of the projection operator P in 
the generator [egPZ + PVA0] indicates that 77,(0 
will be diagonal in /t as required. It is interesting 
to note that, while V e {t) commutes with Z (they 
are functions of atomic position coordinates), the 
operators PZ and PV ( >(t) will not commute in general. 
The formal solution of eq (11) will therefore be a 
time-ordered exponential similar to eq (9). We will 
define an operator 3£ {%) by 



2(g) = eVPZIh 



(12) 



and T' a (t) will be written in an interaction representa- 
tion defined by 



where 
and 



T' n (t) = exp(-it&) U (t (t) 

ih^U a {t) = V{t)U a {t) 

at 



(13) 
(14) 



V(t) = exp(it&)PV e (t)exp(-it&). (15) 
With these new definitions, T„(t) may be written 

T„{t ) = exp (- itHJh) exp (- it g)U«(t), (16) 
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where U„(t) is the (time ordered) solution of eq (14): 

U a (t) = 0expl-{ilh) jV(t')di'Y < 17 > 

2.3. The Line Shape for the Lyman Series 

The line shape function defined by eqs (1-7) and 
(I — 1 1 ) involves a trace over atomic states hence this 
function will describe all spectral radiation emitted 
by the atom. Since we wish to study only one line 
(or a group of overlapping lines) at a time, we may 
pick out the relevant terms in this trace as discussed 
in section 7.2 of I (see eq (1-51)). This selection 
may be accomplished by restricting the atomic dipole 
operator d to have matrix elements only between 
the initial and final states which contribute to the 
line, or lines, of interest (cf. eq (12) of ref. [9]). 

For the Lyman lines in hydrogen there is only one 
final state (neglecting spin), which we have denoted 
by |1); the initial states for any given line are the 
H a eigenstates |/v/) which correspond to the given 
initial energy level E„. All matrix elements of R 
vanish in the ground state, (1|R|1) = 0, hence we 
have (l|f/ a |l) = l (this is why we have chosen the 
Lyman series as an example). 

For convenience we will now specify the states 
\nJz) to be the parabolic states which diagonalize 
both [10, 11] H a and PZ. The diagonal elements of Z 
will be denoted by Z,J and the diagonal elements 
of H a (the energy eigenvalues) have already been 
denoted by E„. There is no loss in generality in choosing 
parabolic states because the line shape is defined 
in terms of a trace over atomic states and we are 
free to choose any complete set of basis states when 
evaluating a trace (i.e., a trace is invariant under 
a unitary transformation of the basis states). The 
atomic density matrix, p" (//„), is also diagonal in 
\/d) and its matrix elements are denoted by p["K 

Using eqs (1-7), (I — 1 1), (1-47) and defining a fre- 
quency variable Aw = (a> — ou„i) for the particular line 
of interest, the line shape function /(co, %) may be 
written in the form 

/(«,*) = (p<?V*r) ReX//'W|d|l> -<l|d|//./> 



£ ; ix, \) = imv 2 + Vf,{x) 



v.v), and V p (x) is the potential of interaction between 
electrons. 

The thermal average is contained entirely in /'(/) 
and the remaining sections are devoted to the evalua- 
tion of this function. When Fit) is known, the line 
shape is readily obtained from eqs (3) and (18). 

3. Further Simplification of the Time 
Development Operator 

3.1. General Purpose 

Although the time development operator was con- 
siderably simplified by the no-quenching approxima- 
tion, the resulting expression, eq (16), is still quite 
complicated and further simplifications are necessary 
in order to evaluate F(t). To facilitate these simpli- 
fications we note that when we use the classical 
trajectories discussed in section 8.4 of I, the inter- 
action potential V ( \t) in eq (2) may be written in the 
form 

MR, x, v, t)=<>R • 2,j9e(x Jn v,, t) =XjVj(t) (22) 



Mx/, v,, t) =e(x j +v j 0/|x i +v j *| 3 ' 



(23) 



where 8V(xj, v y , t) represents the electric field at the 
atom produced by the yth perturber; Vp) is simply a 
shorthand notation for eR -£v(xj, Vj, /), the interaction 
between the atom and the yth perturber. The form of 
the potential, XjVjit), suggests that it may be possible 
to express U (l (t) in terms of a product of time develop- 
ment operators Uj(t) for the individual electron-atom 
collisions. In section 3.4 it will be shown that such a 
product form is obtained without approximation if 
the time of interest is very short. For most cases how- 
ever, the interaction V{t) contains terms like 

Vj(t) = exp (it&)PVj(t) exp (-it ) ) 

which do not in general commute with one another, 
and this prevents us from obtaining the product form 
for U a (t). Nevertheless, if the Vj(t)do not overlap in 
time (or if this overlap is negligible) the ordered 
exponential in eq (17) may be expressed as an ordered 
product: 



[\^\it(^-(^Z^Ifi)]{ /// .\F{t)\ /f /')dt (18) ^ exp 

where 

F(t)=Q- ] f f exp[-£„(x, v)// ( 71c;,(R,x,v,/Wxdv. 



-(M) 



j^V(t')dt'] 



(19) =0nje*p\-(Uh)( VM')dt' 

(20) where Vj{t) is defined by 



(21) 



Vj(t) = exp (it^)PVj(t) exp {-it&). 



(24) 



(25) 



Q= \ \ exp[-£/,(x, v)/AT]dxdv; 

(Equation (24) may be verified by comparing terms in 
E () represents the energy of the Y electrons which are the series expansions.) This result will be the basis 
described by x^ (xi. X2, . . ., x.v) and v=(vi, v-2, . . ., of the following approximations. 
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3.2. The Impact Approximation 

The purpose of the impact approximation is to dis- 
entangle the overlapping collisions [12, 13] in U a (t) 
and permit the use of the product form stated in eq (24). 

To do this, the impact approximation assumes that, 
for collisions which overlap in time, the average colli- 
sion is weak enough that it may be described by a 
perturbation series in powers of Vj(t). Collisions for 
which (Vrlfi) < 1 may be treated this way (cf. sec. 7 
of I). Strong collisions, (Vr/h) > 1, for which such an 
expansion is not practical, are assumed to be well 
separated in time so that two strong collisions never 
occur simultaneously. If a weak collision overlaps 
with a strong one, the weak collision is neglected and 
strong collisions are thus disentangled. The entangle- 
ment between overlapping weak collisions first 
appears in the second order terms in the series ex- 
pansion of eq (24). In footnote 7 of ref. [13] it is stated 
the entanglement in the second order terms vanishes 
because the terms involving different electrons will 
average to zero (recall that U a (t) appears only under 
an average as stated in eq (19)). This result was ob- 
tained in ref. [13] by treating the electrons as statisti- 
cally independent particles (using Debye shielded 
fields to approximate their correlations). Although 
we have not made such an assumption, we note from 
eq (22) that the second order terms involve an electric 
field autocorrelation function (g<>(t)gie(t') ) ; and the 
use of statistically independent quasi-particles is 
known to be a good approximation in evaluating such 
functions. We may thus regard the impact approxi- 
mation as being valid to second order for weak 
collisions. 

Using eqs (17) and (24), we see that, in the impact 
approximation, U a (t) may be written in the form (cf. 
p. 497 of ref. [12]) 



Ua(t) = ^njUj(t) 



Uj(t) = ^exp -{UK) 



\[ww\ 



(26) 



(27) 



where €o is a time ordering operator which keeps 
the collisions in chronological order; that is, the 
Uj(t) which describes the first collision must be the 
first operator from the right in the product YljUj(t) , 
and the second collision must be the second from the 
right, etc. To state the action of c7 more explicitly, 
we note that the variables Xj and Vj in Vj(t) may be 
transformed into the collision variables po, Vo, £o, etc., 
(see appendix) where Jo is a reference time for the 
collision {to is the time of closest approach if the 
collision is completed). The operator &§ simply re- 
quires that to for the first collision always be less than 
to for the second collision, etc. (this ordering is neces- 
sary because the Vj(t) for different collisions do not 
in general commute). It is important to note that, when 
the impact approximation is valid, Uj(t) refers to a 
single collision which may be either weak or strong. 



3.3. The Line Wings 

In the wings of a line, the time of interest 1/Aa> (cf. 
eq (1-52)) is relatively short and it may be possible 
to neglect the exponentials, exp (itST), in V{t). To show 
this, we write eq (14) in the form 



ihj- UX-\U (l (t)\„A-') = ^ exp[itemZ„/-Z„/»)lh] 



dt 



(^\VAt)\^ J, )(^ ,f \Ua(t)\^') 



(28) 



(recall that U a {t) is diagonal in^z). If the time of interest, 
(l/Ao>), is small enough that 



ba>>eg ay (Z„/-Z„/»)h 



(29) 



for all Jy which pertain to the initial level, then the 
exponential in eq (28) may be neglected for£f<^? av , 
where ^av denotes the average ion field strength. For 
?>^ lv , the electron broadening is negligible (as we 
will discuss in sec. 7), thus, if eq (29) is satisfied, we 
may neglect the exponentials involving Z. In this case, 
the solution of eq (28) is 



U„(t)=rexp -(»/*) PV,(t')dt' 



r 



(30) 



and the time ordering is still necessary because PV v {t) 
will not in general commute with PV ( \t'). 

Using the impact approximation, U a (t) is again given 
by the product form stated in eq (26) and Uj{t) is given 
by 



Uj(t)=0exp\-m 



f^pvjb'w}- 



(31) 



Again the time ordering appears because PVj(t) does 
not commute with PVj{t r ). 

For weak collisions the exponential is valid only to 
second order in PVj(t) (because of the impact approxi- 
mation) and it can be shown that, to this order, the 
operator is unnecessary. To show this we note from 
eq (22) that PVj(t) may be written in the form 

PVj(t) = e[PX%At) + PYVjit) + PZ9Jit)l (32) 

It is the orthogonal components of /^R which do not 
commute; that is, PXg \ r {t) commutes with PX%At') 
but not with PY&Jit') or PZ% z (t'). To second order the 
terms which do not commute will contain factors like 
(&At)&y{t r )) (recall that U ( ,(t) appears only under an 
average). For a spherically symmetric (about the atom) 
distribution of electrons these terms will average to 
zero. We may therefore drop the time ordering and 
write 



Uj(t) = exp \-(ilh)i' PVp'W 



(33) 



which is valid to second order in PV,{t) for the line 
wings (as specified by eq (29)). 
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The above argument is not valid for strong collision 
because the series expansion of Uj is not applicable 
We will nonetheless assume that eq (33) provides \ 
good description of the influence of strong collisions 
this approximation will be discussed in section 7. 

3.4. The Quasi-Static Region 



In the far line wings where the time 
1/Acd, is much 
collision, 



less than the duration time 
rAw > 1, 



interest, 

r for a 

(34) 



an interaction 
as dynamic (e.< 



which would normally be regarded 
an electron-atom interaction) will 
take on a quasi-static aspect because the perturber 
motion during the time of interest is negligible. 
The transition from dynamic behavior in the line 
center to quasi-static behavior in the far wings is 
not properly described by current line shape theories 
because these theories approximate U„(t) by the 
first few terms in its series expansion and, as will 
be shown in section 7, this approximation breaks 
down in the quasi-static wings. Nevertheless, it can 
be shown that U„(t) takes a relatively simple form 
in the quasi-static wings. In fact, the product form 
WjUj(t) may be obtained without approximation 
and no time ordering is necessary in the operators 
Uj(t). To show this, we note that the time ordered 
exponential in eq (30) (which is valid in the wings) 
may be expressed as an exponential involving a series 
of commutators [14]: 



(35) 



U ft (t) = exp 

+tmv 



dlfi) PV v {t')dt' 



For quasi-static interactions we have V—fi\io (cf. p. 
23 of ref. [4]) and t 1 P V ;,{{)) I bfi ^ 1/tAw<^ 1. Thus, if 
eq (34) is satisfied, the commutator in eq (35) is neg- 
ligible and U a (t) may be expressed in the product 
form as stated. 

4. Series Expansions of F(t) 

In this section we will obtain series expansions for 
the function F{t) defined by eq. (19). The various 
methods of handling these expansions give rise to 
different types of theories (e.g., the one-electron 
theory, impact theory, etc.). 

The classical Hamiltonian /^(x, v) appearing in 
F(t) was given, in eq (20), by 

£ p (x, x)=^mvjl2+V p (x u . . .,x.v) (36) 

3 

where Xj and vj refer to the ;th electron and V p (x) 
represents the electron-electron interaction. We may 
define position and velocity distribution functions 
P(x) and W(v) by 

P(x) = exp {- V p (x)lkT}l J dx exp {- V p (*)lkT) (37) 

W(\j) = exp {- mvjl2kT}/ J d\j exp {- mvj/2kT} (38) 

W(v) = YljW(vj), (39) 
and these have the property that 

r(v)P(x) =Q~* exp {-£ p (x, v)/kT} (40) 



where Q was defined by eq (21). Using the product 

\PVAt') PVAt")](lt' dt"+ \- f<>rm f ° r U(l ^ ^ iven in eq (26 ^' We ()l)tain 

F{t) = aAdxP{x) llljdvjWivjWjit). (41) 
If the commutator is negligible, PV v {t) may be 

regarded as an operator which commutes with itself Following the technique used by Baranger and Mozer 

at different times; similarly we may regard PV,{t) [15], we define a function <p(R, xj, t) (which should 

and PVj(t') as commuting operators for all i 9 j\ t< not be confused with the wave functions if used in 1) by 
and t' . In this case, U a {t) may be expressed in the r 

product form ILjUj(t) where Uj(t) is given by eq (33). 1 + ?(R, *j, = I dvjW(vj)Uj(R 9 x,-, Vj, t) . (42) 

To show when the commutator is negligible we 

expand V e (t) about / = (in the far wings t is a small Using <p(j) as a shorthand notation for p(R, xj, t) we 

parameter). To lowest order the integral over PV ( >(t') have 



is given by t PV e (0)lh and the integral over [PV e (t') % 
PV r (t")] becomes f" [PF,(0), PV' t .(0)]l6h\ where 
V',.(0) denotes (dV e (t)ldt) at / = 0. Since the magni- 
tude of [PV f .(0), PV' f .(Q)] is less than or equal to 
PF r .(0)PK^(0), we need only compare the magnitude 
of (t-PV' ( ,(0)l61i) with unity. To estimate this quantity 
we use eqs (22) and (23) and we obtain 

(BVel<H)=(dx/dt) (dVelfix) - Vo(dVeldx) - (VoVJfh)), 

where p<» and t/« denote the position and velocity at 
the time of closest approach (when PV[ is largest). 
Usingr= (po/f<>) and t— ( 1/Aco). we obtain 

t 2 PV;,(0)/6h - (V r lhAco)(llTAoj). 



J j<k 



(43) 



Substituting eqs (42) and (43) into (41) we obtain the 
series 



F(t)- 



(44) 



^*to =£*<*< <} J<p(jt) • • ■ ^(/if)P(x)rfx. (45) 



where 
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Since P(xi, . . ., x.\) is invariant under an interchange 
of particle coordinates we have 

F XI (t) = [N\IMl(N-M)l] 

J . . . J <p(l) . . . <p(M)P(x, x. v )dx, . . . dx A 

= [N\IMl(N-M)l] 

. . . \ ip{V) . . . <p(M)/Vxi, . . ., x M )dxi . . . C?XM, 

(46) 

where the reduced distribution function Pa/(xi, . . ., 
xj is defined by 

P M (xi, . . .,x M )=[. • J P(x 1? . . ., 

xm, Xiif+i, . . .,x A ')dx.u + i . . . rfx.v. (47) 

Pm(xi, . . ., x.\/)<ixi . . . Jx.y represents the proba- 
bility of finding the M particles labeled 1 through M 
in the volumes dxi, . . . rfx.u at the points xi, . . .,x.u 
(for a further discussion see ref. [16]). The com- 
binatorial factor in eq (46) represents the number of 
ways for choosing M electrons from the perturbing gas 
of N identical electrons. 

The reduced distribution functions P.u(xi, . . ., xa/) 
may be expanded in an Ursell expansion [17], the first 
three terms of which are given by 

rp 1 (x 1 )=gi(x 1 ), 

r 2 P 2 (x 1 ,x 2 )=g 1 (x 1 )g 1 (x 2 )+g- 2 (x 1 ,x 2 ), 

r 3 P 3 (x 1 ,X 2 ,X3)=^l(x,)g-,(x 2 )g-i(X3) 
+ gl(Xi)g 2 (x 2 , X3)+#l(x 2 )g 2 (Xi, X;0 

-h^(x,)^(x,, x 2 ) + &(xi, x 2 , xO, (48) 

where T is the volume of the system; the general term 
in the Ursell expansion is given by eq (9) of ref. [17]. 
The function gi/^is the one-body probability function; 
g,(xi)^i(x 2 )/F 2 would be the two-body probability 
function if there were no correlation between particles. 
Since there is a correlation between the particles at 
xj and x 2 , this correlation is represented by the two- 
body correlation function g" 2 (xi, x 2 ). Similarly P :i is 
expressed by a group of terms which contain no more 
than two-body correlations plus gbCxi, x 2 , x.O, the 
three-body correlation function. The general term in 
the Ursell expansion expresses Pa/(xi, . . ., Xm) as 
a sum of terms which contain no more than (M-l)-body 
correlations plus £a/(xi, . . ., x.u), the M-body correla- 
tion function. Following Baranger and Mozer [15] we 
define a function 



and, after considerable algebra, it is found that F{t), 
as given by eqs (41) and (43), can be expanded in terms 
of the gM and then resummed to yield 



F(0 = ^<>exp \ 2 (n M IM\)h M (t) 



(50) 



h M (t)=j <p{l) • • • <p(M)g M (xu • ■ ., 

Xa/)g?Xi . . . G?X\/, 



(49) 



In the derivation of eq (50) it was necessary to let N 
and V become infinite in such a way that the density 
n = NI'V remained constant. In this new series for 
F(t), each successive term represents a higher order 
correlation effect. 

We now have two entirely different expressions 
(eqs (44) and (50)) for F(t). To determine the utility of 
these expressions, we will briefly outline the physical 
significance of the terms in the series expansions in- 
volved. To understand the physical significance of 
F(t) itself, we note that, except for multiplicative 
constants (e.g., d matrix elements), F(t) is the same as 
the autocorrelation function C(t) discussed in section 
7 of I. 

The Fj{t) term in eq (44) describes the electron- 
atom interactions (collisions) as though they are com- 
pletely independent of one another. From eqs. (42) 
and (46) we see that F^t) is given by a sum of terms, 
each of which contains only one Uj operator, thus, to 
this order the effects of the interactions are simply 
additive. The terms F 2 , . . ., F\ account for the fact 
that the effect on F(t) of one interaction may be in- 
fluenced by other interactions. When this influence is 
taken into account, the effects of the collisions are no 
longer simply additive. This influence of one inter- 
action on another comes about in basically two ways. 
The most obvious influence is due to the correlation 
between the perturbing electrons which are de- 
scribed by the correlation functions g\i for M^2. How- 
ever, even when these correlations are ignored (by 
setting gM = for M ^ 2), the F 2 , . . ., Fm terms do 
not vanish and it is still possible for the outcome of 
one collision to be influenced by other collisions. This 
remaining influence comes about when a collision, or 
a series of collisions, produces a large change in the 
state of the radiator thereby reducing the autocorrela- 
tion function F(t). If the state of the radiator is 
thus materially altered from its original value, the 
effect of subsequent collisions will be correspondingly 
modified. For example, there is a high probability that 
a strong collision will cause an inelastic transition 
between states having different principal quantum 
numbers and it is obvious that such a transition cannot 
be neglected when calculating the effect of the next 
collision. Since it takes several weak collisions to 
build up such an influence, we expect that F \{t) will 
provide a good approximation to F (t) for times which 
are shorter than the mean free time between strong 
collisions (for those cases where a strong collision 
does occur during this time, the effects of subsequent 
weak collisions are negligible in comparison). For 
longer times of interest it is necessary to consider all 
of the terms Fi, . . ., F\, and in this case eq (50) is 
more appropriate than eq (44). 
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In eq (50), the hj(t) term describes the electron- 
atom interactions as though the position of an electron 
is not influenced by the presence of the other elec- 
trons. That is, the electron-electron correlations are 
neglected in this term (although the remaining influ- 
ence of one collision on the effects of subsequent col- 
lisions is properly included). The h 2 term accounts 
for pair correlations between electrons, and subse- 
quent terms describe higher order correlations. Since 
it is usually not necessary to consider more than pair 
correlations, it would seem that eq (50) is much more 
useful than eq (44). This would indeed be the case if 
we could always make calculations based on eq (50); 
as we will show in the following sections, the calcula- 
tions based on eq (50) become quite difficult for short 
times of interest and in this case eq (44) is of more 
practical value. 

As a further comparison of eq (44) and (50), it is 
interesting to note that, in the microfield theory 
[15], the approximation F — exp {nh\) gives rise to the 
Holtsmark distribution function (no correlation ap- 
proximation) whereas the approximation F — \-\-F\ 
gives the asymptotic Holtsmark function (or nearest 
neighbor approximation). This indicates that the 
approximation F(t) — ^o[l +F\(t)~\ is just the small 
t limit oCF(t)^ffo exp {nhi\t)}. To verify that 
this is the case, we note from eqs (42) and (49) that as 
t— >0, h\(t) — > and exp {nh\(t)} — * 1 + nhi(t); sub- 
stituting the functions [16] Pi(x\) =njN and^i(^i) = 1 
into eqs (45) and (49) we see that F\{t) =nh\(t), thus 
\+F\(t) is indeed the small t limit of exp {nh\(t)} . 

5. The One-Electron Theory 

The one-electron theory is designed to provide 
a description of the wings of a line profile. The version 
of this theory presented in this section will be shown 
(in sec. 7) to provide a consistent description of 
a line profile from the half width to the quasi-static 
wings. This theory is based on the approximation 
F(t) =*&o[l + Fi(t)] which is valid for small times 
of interest or large Ao>. Since Fi(t) is given by a sum 
of terms each of which contains a single Uj(t) operator 
(rather than products of the Uj), the chronological 
ordering operator Gq will have no effect on F\. Using 
[16] Px(x { )=nlN m eq (45), we obtain 



■/ 



F(t) = l + n <p(R,x,, /)</x, 



(51) 



Notice that this expression involves the average 
of the time development operator for a single colli- 
sion; this is the origin of the name "one-electron" 
theory. 



The first term in F(t) (i.e., the constant 1) will 
give rise to a delta function, 8(Ao> — egfZ,,///*), when 
it is substituted into eq (18). Performing the ion micro- 
field average, as required by eq (3), it is clear that 
the contribution to J> (co) resulting from this term 
will be proportional to & [hAcoleZ,,./). This represents 
the line shape which would be produced by the ions 
alone. 

The influence of the electrons (as well as some 
electron-ion coupling) is contained in the second 
term in eq (51): 

Fi(t) = nl<p(R, xi, t)dxi. (52) 

Using eqs (33) (which is valid only in the line wings) 
and (42), and noting that I W{v\)d\\ = 1, we obtain 

F,{t) = iAdx x W(i h ) [r/x, 



f 



exp\-(i/h) PVi(R,in + Yit')dt 



1 



(53) 



As discussed in the appendix, we may transform 
to the more familiar collision variables (p , vo, to) 
where p and vo denote the impact parameter and 
electron velocity while to is some reference time 
in the collision (the time of closest approach in a 
completed collision). In terms of these variables 
we have 

F 1 (t) = nld€ll vW(vo) dv Q I podp I dt 

exp {- (ilh)j PVx (R, po, vo, t' + to)dt'\- 1 (54) 

where 

Fi(R,po,vo, t' + t )=e 2 R 

• |>o + vo(*' + * )]/[p8 + t>8(*' + to) 2 ] 3/2 . (55) 

The Euler angles 11 denote the orientation of the 
"collision axes" relative to the direction of R, the 
atomic dipole. We may therefore interpret the integral 
over fl as an average over all possible orientations 
of the electron trajectory [5] relative to the direction 
of R. 

To make full use of the collision axes, we take the 
matrix elements of Fi(t) as required by eq (18). 
These matrix elements are given by 



U\F X (t) \n£) = n j dil P vNTMdvo P podpo f * dt {) 



,/■ 



exp l-(ilh) JPF,(R, po, vo, t' + t {) )dt'\-\ 
411 



/')• (56) 



We may now perform a rotation of the "atomic axes," 
through the Euler angles H, so that R will point in the 
same direction as po, and v () defines the x direction. 
The rotated wave functions will be given by [18] 
^(n)|/ii/), where 2 (O) is a rotation operator. The 
interaction potential referred to collision axes is 



= e 2 [Zp () +^o(^ + /o)]/[p(H^(/' + ro) 2 ]^, (57) 

where X and Z denote the x and z components of R, 
and we have the identity 



iJk\ 



exp{-(i/ft) jPVi(R 9 po,vo 9 t' + to)dt'\-l M 



= LnA\9- l (0) 



exp \-(ilfi) J PF r (Z,Z,p , i*>, *' + to)*')-l ^(ft)M'V (58) 



Substituting this identity into eq (56), we obtain 
where 



(59) 



FV(0 



/? I t$JF(i;o)di;« I podpo cfeo 



exp ]-(£/*) | PFc(*,Z, p«, i*, *' + *,)&' f-1 

(60) 



The average over (1 in eq (59) may be readily performed 
by transforming to spherical wave functions [19] 
and using the properties of .(/ (£1) as given by Ed- 
monds [18]. The function F\{t) may be evaluated by 
numerical methods, and the line shape is then obtained 
by substituting the result into eq (18) and performing 
the transform. Further details of such a calculation 
together with numerical results will be reported in a 
future paper. 



6. The Impact Theory 

To obtain the results of the impact theories we use 
the Uj(t) which were obtained by the impact approxi- 
mation, eq (27), and the approximate calculation of 
F(t) is based on its expansion in terms of correlation 
functions, eq (50). For most plasma broadening prob- 
lems one is interested in a temperature-density range 
where three-body correlations are negligible, hence 
it is necessary to consider only h\(t) (no correlations) 
and h 2 {t) (two-body correlations). In the impact 
theories, electron correlations are not treated with 
correlation functions (such as ^(xi, x 2 ) which ap- 
pears in hi(t)) . Instead, in these theories the 
electron correlations are approximated by impact 
parameter cutoffs. Shielded field arguments are 
used to derive a cutoff of approximately A./>, the 
Debye length, in the range of the electron-atom inter- 
action [20|, and the frequency dependent Lewis 
cutoff [21] (i; av /Ao)) is used in place of Xn when 



(t; av /AG>) < A/; (this will be discussed in sec. 7). A 
similar cutoff behavior has also been obtained, by 
means of correlation functions, within the framework 
of a quantum mechanical theory [22]. However a 
more accurate treatment of the electron correlations 
[23] has shown that, while the cutoff at (v av /Aa)) is 
valid for Aoo > co,,, the cutoffs in the line center are in 
error. The cutoff at the Debye length (or at 1.1 Kd in 
some theories) should be multiplied by (2.718)" 1 - — 0.6 
resulting in a change on the order of 20 percent in the 
broadening function <$>„!, used by Griem et al. 
[2, 9, 13, 20]. For Aoj — oo/„ there is a possibility of 
further corrections due to dynamic correlation effects. 
To derive the results obtained by the impact theory, 
we will assume that the effect of electron correlations 
may be represented by some correctly chosen impact 
parameter cutoffs. We thus need to consider only the 
function h\(t) as given by eq (49): 



F(*)~0oexp{iiMO}- 



(61) 



For a fluid [16] such as a plasma, g\ (xi ) = 1 and h\ is 
given by 



Ai (0 =jj rfxidvir(v,)[i/,(R,x lf v,, 0-1]. 



(62) 



We note in passing that if the exponential in eq (61) 
is expanded, the first two terms in the series are iden- 
tical with the approximate expression for F(t) , eq (51), 
which is used in the one-electron theory. A further 
discussion of this correspondence is given in the 
following sections. 

Using eqs (25) and (27), h\ (t) is given by 

M0=[ fdK«hiW(vi)[0exj>{-(ilh)( e"'* 

PF,(R,xi,v,,Oe-' rir <fc'}--l]- (63) 

Transforming to collision variables (see appendix) 
and denoting the average over p , ^o, and D by a sub- 
script av, this becomes 
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h 1 {t)=\ X dt {) ^exp|-(/M) j* e lt '* PVi(t' + t Q )e- it '*dA-~l . 



Changing variables from t' to s = t' + tQ and using the 
identity exp [e~ A B A ] = e~ A e B e A , we obtain 



hi(t)=j cfloexp(-fto^)r^exp|-(£M) j " e^PV l (s)e-^ds\-l 



exp (ito3Z). 



(64) 



(65) 



We assume that a collision time t exists such that 
Vi(s) = unless ^ s ^ r. The completed collision 
assumption, which is an important part of most 
impact theories, states that if a collision occurs, 
it can be completed. That is, if some value of s in 
the interval [to, to + t] satisfies O^s^r, then the 
completed collision assumption states that the col- 
lision should be completed in that interval, i.e., 
r<>^ 0^ s.^t^ (*o + f). Using this restriction on h 
and extending the limits U<>, / + /<>) on the s integral to 
(— oo, + ») (since V\(s) =0 outside the range ^ s ^ 7 
anyway), eq (65) becomes 



h } (t) = dt {) exp (— ftoSf) [S— l] av e *p (iuS) 

J-U-T) 

cu-t) 

= dto exp (it(&)[S— l] av ex P (— i£o<3T), 



(66) 



where S is an .S-matrix for a single collision (weak or 
strong) defined by 



I 



S = ^exp \-(ilh) e»*PV x (s)e " J ds 



(67) 



The completed collision assumption will provide a 
good approximation to h\(t) since r is usually very 
small (i.e., t > r), however, if t < r, the collision cannot 
be completed and this approximation breaks down. 
Using the time of interest 1/Aco from eq (1-52), we may 
state that the completed collision assumption will be 
valid when 

tAco<^1, (68) 

where r is some "representative" collision time [3]. 
This inequality is sometimes used as the validity cri- 
terion for the impact approximation [3]: this point 
will be discussed further in section 7. 

Substituting eq (66) in eq (61) and using the operator 
identity 



o exp 



i exp 



-t,)A)B exp (uA)dt» 



(69) 



we obtain 



F(*)=e«^)^exp [-i(t-T)&+(t-T)®< t ] 

= e ««-*>* exp [+JTJ-A] exp [-i^ + f4> M ], 



where (cf. eq (42) of ref. [3] or eq (28) 

<D„ = n[S-l] av . 

ise 

tO„ < 1 



II 



which is simply a statement that the impac 
tion is valid (as we will show in sec. 7), 
becomes (ci\ eq (39) in ref. [3] or eq (30) 

F(t)=e i <* exp [-i^ + *<&,,]. 



of ref. [12]) 

(71) 

(72) 

t approxima- 
tion eq (70) 
of ref. [12]) 



(73) 



Substituting this F(t) into eq (18), we obtain 
I(oj^) = (p { :: ) Itt) Re 2 <«*|d|l) -(l\d\nA') 



//■' 



f 



M| exp [i(ba)-egPZIh)t + t&u]\»*')dt, (74) 



(70) 



(Comparing with eqs (9) and (32) of ref. [13], we see 
that this result is formally identical with the results 
obtained by the impact theory for the Lyman series. 

It is important to note that, in most calculations of 
hydrogen line broadening, the exponentials exp 
(is3?) in the integrand of the S-matrix (eq 67)) are 
either ignored or approximated by a cutoff [24|. 
To show how this cutoff is obtained, we note that when 
e&r(Z„/— Z„/>) lh < 1 the exponentials may be re- 
placed by unity. When e&T(Z„/—Z„/')lh^>l n the 
exponentials oscillate rapidly and there will be only a 
negligible contribution from this region of the integral. 
Defining a frequency 

o>«M „/> = e9(Zu-Z,J')lh (75) 

and noting that r = po/zo, we define an impact param- 
eter cutoff Pmatf^tWct),,^/,/' (in ref.[24]co // / /// /' is called 
Aco.s). When this cutoff is used, the exponentials exp 
(is 3£) in eq (67) are replaced by unity; this cutoff 
reduces the integrand to zero when rco,,/ \ „/,• > 1 , 
thus approximating the effect of the rapidly oscillat- 
ing exponentials, and when toW,/,/'< 1 the exponen- 
tials are effectively replaced by unity. This cutoff 
is used simply as a method of simplifying the numerical 
calculations and it could be corrected without any 
modification of the theory itself. 



413 



7. Discussion of Approximations 

7.1. General 

In the preceding sections we have derived the 
familiar results of the impact and one-electron theories 
using the classical path formalism developed in I. 
The advantage of these theories lies in their simple 
classical description of the perturber dynamics; their 
disadvantage with respect to more recent quantum 
theories [6, 22, 25 ] , is in the large number of approxima- 
tions which are necessary to achieve this classical 
description. Since an understanding of these approxi- 
mations plays an important role in the application and 
improvement of classical path theories, these approxi- 
mations will be discussed in some detail in this 
section. 

7.2. Preliminary Approximations 

In our treatment of both the impact and the one- 
electron theories we have used a dipole approximation 
to the perturbation potential, a static treatment of the 
ion field, and the no-quenching approximation. 
It should also be mentioned that we have not con- 
sidered any perturbation of the final state in the radia- 
tive transition because there is no linear Stark effect 
in the ground state of hydrogen and any higher order 
ground state perturbations (e.g., quadratic Stark effect) 
are unimportant [26|. For lines other than those in the 
Lyman series, lower state perturbations may have to 
be considered. 

To show that the static approximation provides a 
good description of the ions, we consider the change 
8& in an electric field &=ejf 1 during a time 67. Using 
8gr^8t(d&ldt)^v8t(d%ldr) - (v&8tlr). we see that 
(S^f/^f) may be estimated by (dt/r) where r is the dura- 
tion of a collision. If 8t is replaced by the time of inter- 
est (1/Aw). we see that ( SgT/ST ) will be negligible when 



r Aw 5> 1 , 



(76) 



as mentioned earlier. If a perturber is quasi-static 
it produces a Stark splitting hAoj given by hAa)=V 
(see p. 23 of ref. [4]). Using V ' — Ze 2 Rolp$ for an ion 
of charge Ze, this means that pi — Ze 2 Ro/hAa>. Since 
the velocity distribution is sharply peaked about 
v av , we estimate t by (po/v ay ) and eq (76) becomes 



Aco> Aw r = (hvlJe 2 R {) ). 



(77) 



The critical frequency Aco r is sometimes referred to as 
the Weisskopf frequency. Since Ao> c for ions lies al- 
most at the center of most lines, the static ion approxi- 
mation is usually justified over most of the line profile 
(however see also ref. [27]). In section 3.E of ref. [6], 
it is shown that the effect of a static ion field may be 
described by the microfield average stated in eq (3). 
In the no-quenching approximation it was assumed 
that the 77,(0 matrix elements between states having 
different principal quantum numbers could be ne- 
glected when calculating C(t). These "off-diagonal" 



matrix elements represent collision induced transitions 
which transfer the excitation of the atom from one 
energy level to another. In section 7.2 of I, it was shown 
that collision induced transitions between excited 
states correspond to switching the mode of oscillation 
in a classical oscillator. It was also shown that this 
mode switching effect could be neglected if the energy 
separation of the excited states is much larger than the 
halfwidth of the lines being studied. In hydrogen, states 
having the same ^are degenerate (or nearly degenerate 
in the presence of an ion field) and collision induced 
transitions between these states are very important; 
this is the reason for regarding all states with the same 
/n as a group of initial states. Low-lying states with 
different <n are well separated in energy (much more 
than the halfwidths of the lines) hence the transitions 
described by {n^\T' a (t) \<nA') rnay usually be ne- 
glected. However, for large <n the levels are closer 
together and some of the higher series members may 
overlap; in this case it may be necessary to include 
states with different <n in the group of initial states. 
In section (2.2) the no-quenching assumption was 
made by setting ( < nJ>\T[ l (t)\'riA!) to zero. This neglects 
the quenching terms which would be added to eq (18) 
as discussed above, however, it also neglects the in- 
fluence which (^iA\T' a {i)\^i!A f ) has on {riA\T' a \nA!). 
That is, even if a line is isolated in the sense that the 
quenching terms which would be added to eq (18) are 
negligible, it may still be necessary to consider the 
influence of the inelastic transitions \nA)^\n'A') 
when calculating {nA\T' ( \nJ:'). To show when this 
influence is also negligible we consider the differential 
equation for T' (l (t), eq (6). This equation may be writ- 
ten in the form 



ih-{nA\T (l {t)\W) 



A" 

+ 22 ex P {it<*)\*s<"} (<n^\e > &Z 



■v e (t)\n'W)ww\n(t)\^ f ), (78) 



where co,,,/ •= {E„ — E„")lh. From this equation we see 
that setting {,J\T ( \„'/f') to zero simply neglects the in- 
fluence which the "off-diagonal" matrix elements of 
[egZ+V ( ,(t)] may have on {„/-\T a \,J' ) . If eq (78) 
is solved in the usual iterative manner [7] by inte- 
grating over t . the term exp (ita)„„") will oscillate rap- 
idly and contribute essentially nothing for t >(l/co /// /). 
This indicates that the (n£\Z\n&) terms will be 
negligible compared with the (,dt\Z\nA') terms un- 
less t < (1/co „,/). Since the time of interest, 1/Aoj, 
is much larger than (l/o),,,,-) for isolated lines, we may 
neglect the terms involving "off-diagonal" Z matrix 
elements. The same results apply for the terms in- 
volving "off-diagonal" elements of V v {t) unless 
t > r. If / > r, a similar argument indicates that the 
"off-diagonal" V e (t) terms are negligible only when 
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tq) /// /> 1. Iii this case we use T — po/v av (the velocity 
distribution is sharply peaked about v in ) and we sec 
that the condition TCt>>,„" > 1 is violated only when 
P» < (hv av l&E„„") where A£7,,/' = {E„— E„"). Substi- 
tuting this value of p {) into (Vr/h) — (e 2 /?oMp«i>av) 
(recall that V — e 2 Ro/pl for a hydrogen plasma) we 
find (VTlh)~{e 2 RiM,,, t "lh 2 v\ v ). Using v 2 v =* (MT/m) 
and /?o =* (3n 2 ao/2), where ao = ft 2 /me 2 , this becomes 
(VTlh)>(ri 2 AE„,rl2kT). If we use A£_ »=£ X |(1M 2 
— (1///) 2 |, where £,. = 13.6 ei;. is the ionization energy 
of hydrogen, we find that the off diagonal-terms are 
important only when 



(Vr/h) >(EJ2kT)\ l-(/v/,/') 2 



(79) 



A Since E x ^> ATone usually has (EJ2kT)\l- U/>/T| >1 
and, in this case, eq (79) indicates that the (///£" 1 77, |///') 
terms in eq (78) are important only during strong 
collisions (in which case T' n averages to zero, as 
discussed in sec. 7 of I, so no error results). For 
some higher series members (very large " and //') 
it can happen that (EJ2kT) |1 - U/**") 2 | < 1, and 
eq (79) can be satisfied even for weak collisions. In 
this case, the "off-diagonal" terms in eq (78) must 
be retained when calculating (/*/|77,|,//'). To see 
when this occurs, we consider the case where 
//' = ,*±1 and „> 1. In this case the quantity 
(EJ2kT)\l-(*l,i") i \ will be less than unity when 



n>(EJkT). 



(80) 



Thus, for lines involving v*> (EJkT). we may not be 
able to neglect the influence of inelastic transitions 
hetween states having different principal quantum 
numbers when calculating ( /t /-\T' a \ ,</■'). This effect 
is usually negligible for hydrogen lines from labora- 
tory plasmas because, for these plasmas, the lines 
generally begin to merge (overlap) when eq (80) is 
satisfied. For some astrophysical applications [27] 
(and for lines other than hydrogen lines) this effect may 
be very important. 

When the differential equation for the matrix ele- 
ments of 77,(0 in the no-quenching approximation, 
eq (10), was written as an operator equation, eq (11), 
it was necessary to introduce a projection operator 
P that picks out the part of an operator which is 
"diagonal" in //. Although P operates on Z and V ( >(t), 
this does not imply that the "off-diagonal" elements 
of these operators are small (indeed they are about 
the same order of magnitude as the "diagonal" 
elements). The appearance of P in eq (11) merely 
indicates that the "off-diagonal" parts of Z and 
V r (t) have a negligible influence on T ( [ in the no- 
quenching approximation. The influence of this 
projection operator is implicit in most line broad- 
ening theories since matrix elements of Z and V are 
usually taken only between states having the same 
principal quantum numbers; we have introduced 
P explicitly in order to avoid errors in handling T,',(t) 
and to clarify subsequent approximations. For ex- 
ample, the only operator in V e {t) is the atomic posi- 
tion operator R and one may therefore commute 



Z with V r (t); this commutation is not possible with 
PZ and PV e (t) , in fact, the matrix elements of the 
commutator [PZ , Pl\.(t)] are about the same order 
of magnitude as the matrix elements of PZ and 
PV ( \t). If the influence of P is only implicit, it is 
possible to make serious errors by commuting Z with 
V ( >(t) when these operators appear in 77,(0 • 

Since PZ does not commute with PV e (t) , the formal 
solution, 77(0 » of ec i (11) is a complicated time- 
ordered exponential. Although this operator has a 
rather simple form for the quasi-static region of the 
line wings (see sec. 3.4), it was necessary to introduce 
simplifying approximations in order to evaluate the 
rest of a line profile. The first of these was the impact 
approximation which assumes that most collisions can 
be described by a second order perturbation expansion 
of Ua(t). Strong collisions (Vr/h) > 1, for which this 
expansion breaks down, are treated by some other 
method such as an impact parameter cutoff. With 
this approximation, it is possible to write 67,(0 i n the 
product form C?oIljUj{t) where each Uj(t) is the time 
development operator for a single electron-atom col- 
lision. It is difficult to give a useful validity criterion for 
the impact approximation; the validity criterion which 
is usually given [3], rAo> < 1, is based on the com- 
pleted collision assumption which often accompanies 
the impact approximation when it is used in the impact 
theory (see sees. 6 and 7.3). Since the impact approxi- 
mation is also used in the one-electron theory, we have 
drawn a distinction between the impact approximation 
and the impact theory. The impact approximation, as 
we have stated it, will be valid when strong collisions 
are well separated in time. If we define a strong col- 
lision frequency, v 8 , and a maximum duration for a 
strong collision, r«, this condition may be given by 



T 8 Vg < 1 . 



(81) 



If all collisions were strong, v s would be on the order of 
the half width, Acoi/2, of the spectral line. For hydrogen 
plasmas we have v s < Acoi/ 2 . The duration of a strong 
collision may be estimated by using (Vr/h) — {e 2 Rol 
hpov ay ) > 1 or p < (e 2 R { )lhv av ) which gives T« =* (polv av ) 
^ (e 2 R { )lhv 2 av ) = (1/Acoc), where Aco r is the Weisskopf 
frequency for electrons. Using v 2 dv — (3kT/m) and 
/?o— (3^ 2 a /2) where ao= (h^rne 2 ) , we have r H 
<{ri 2 h\kT) and T s v s <// 2 {hA(D\/2lkT). Since ^i 2 (hAa>i 12/ 
kT) < 1 for most hydrogen plasmas, we expect the 
impact approximation to be valid for the electrons. 

7.3. Discussion of the Impact Theory 

In the usual derivation of the impact theory [2, 3] 
one calculates the change in (U a (t\0)) during some 
time Aj (recall that U a (t) = U a (t\0)). This change is 
given by 

MU a (t\0)) = (Ua(t + At\0)-Ua(t\0)) 

= ([U a (t+At\t)-l]Ua(t\0)). (82) 

It is assumed that: (1) A£ is large enough that the two 
factors on the right are statistically independent and 
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may thus be averaged separately, and (2) At is small 
enough that the difference equation, eq (82), may be 
replaced by the differential equation (see eq (4.30) 
of ref. [2], eq (38) of ref. [3] or eq (2.8) of ref. [20]) 



dt 



(U a (t)) = [exp (itHalh)<S> (l exp (-itH (l /h)](U«(t)) 



(83) 



which is then solved for (U„(t) ) (<& n is given in eq (71)). 
In this derivation, assumption (1) requires At > r and 
assumption (2) requires At < t (as well as the average 
collision being weak; cf. pp. 508 and 529 of ref. [3]); 
together these assumptions require t > r which is a 
mathematical statement of the completed collision 
assumption (cf. eq (68)). Thus we see that the com- 
pleted collision assumption is an integral part of the 
usual derivation of the impact theory and the removal 
of this approximation would require a complete revision 
of the derivation. This is the reason one frequently 
sees the condition t^> r or tAco <^ 1 given as the validity 
criterion for the impact approximation [3]. 

In our derivation of the impact theory we first used 
the impact approximation to write U a (t) in the product 
form ^olljUj(t) and we obtained a consistent correla- 
tion function expansion, eq (50), for (U a (t)) =F(t). 
To obtain the usual results of the impact theory, we 
assumed that the electron-electron correlations could 
be approximated by impact parameter cutoffs, we used 
the completed collision assumption and we assumed 
that r4>a < 1. It was also noted that numerical calcula- 
tions are frequently simplified by approximating the 
influence of the ion field splitting on the electron 
broadening (i.e., the factors exp (is££) in eq (67)) with 
an impact parameter cutoff p max =volai,J,,W which 
depends on the ion field strength. While this approxi- 
mation does not affect the region of validity of the 
theory, it would nonetheless be desirable to improve 
on this approximation to insure more accurate calcula- 
tions. It is also interesting to note that the condition 
rOr/ < 1 is simply a manifestation of the requirement 
that the impact approximation is satisfied (this was 
discussed in sec. 7.2). To show this, we use eq (71) and 
estimate the magnitude of O r/ by 

<&„ = n[S -l] av - n(Vr/h) av 

where the subscript av refers to an integral over the 
collision variables p , fo, and 11 (cf. eqs (64) and (A. 5)). 
Since 27rnpodpi)v'f ) W(vo)dvodCl = dv where v denotes the 
frequency of collisions described by p^ vq, and fl, we 
see that n(Vrlh) av = J (Vr/h)dv= (V/h) where V 
denotes the average value of V. The statement t<I>« < 1 
now becomes (Vr/h) < 1, which is satisfied if the 
average collision is weak and if strong collisions 
(Vrlh > 1) are well separated in time. This shows that 
the condition T<f> (l < 1 is always satisfied when the 
impact approximation is valid. 

One potentially useful feature of our approach is 
that it permits a consistent improvement of the impact 
theory within the framework of the classical path 
methods. That is, the impact parameter cutoffs, the 



completed collision assumption and the second order 
expansion of Uj(t) were needed only to simplify the 
calculation of eq (50) and it would be possible to cor- 
rect these approximations without having to rewrite 
the theory. In fact, the impact parameter cutoffs have 
already been corrected to some extent by means of 
a more consistent treatment of correlation functions 
[22, 23]. 

The correction of the completed collision assump- 
tion is a matter of some importance because, if we 
estimate the collision duration time by (\z>/i; av )=(l/a)p), 
the region of validity for this approximation, eq (66), 
is given by 

Ao) < o) p . (84) 

That is, an impact theory with a completed collision 
assumption can be used only for values of Aco smaller 
than the electron plasma frequency. It has been shown 
that the completed collision assumption may be cor- 
rected, in the line wings, by means of the Lewis cut- 
off [21]. This cutoff may be obtained from eqs (61) and 
(62) by expanding Fit) to second order in V(t) and 
performing the Fourier transform as stated in eq (18) 
(cf. eqs (6) and (17) of ref. [21]). 

Since the derivation of the Lewis cutoff is based on 
an expansion of F(t) which is valid only in the line 
wings, the use of this cutoff in the impact theory raises 
some doubts as to the validity of such a "modified" 
impact theory in the line center. These doubts may be 
dispelled to some extent by means of a comparison 
with the results of the quantum mechanical relaxation 
theory which does not make a completed collision as- 
sumption. In reference [22] it is shown that the 
results of the relaxation theory essentially reproduce 
the behavior of a Lewis cutoff, and these results are 
not restricted to the line wings. In spite of this useful 
comparison, it would nonetheless be desirable to con- 
firm this behavior within the framework of a classical 
path theory. In this regard, it is important to note that 
the completed collision assumption is not an inherent 
part of our derivation of the impact theory. This ap- 
proximation was introduced simply to compare with 
the results of the usual derivation. It would thus be 
possible, in principle, to calculate F(t) from eqs (61) 
and (62) without making a completed collision assump- 
tion. This has in fact been done, to second order, 
[28], and it is found that the results are identical with 
the results of the relaxation theory. It thus appears 
that a "modified" impact theory which does not make 
a completed collision assumption (such as ref. [28]), 
or a theory which corrects this assumption by means of 
a Lewis cutoff, will have a wider range of applicability 
than that stated in eq (84). 

If the completed collision assumption is not made 
(or if it can be corrected by some approximation such 
as the Lewis cutoff), the only serious limitation on 
the applicability of such a "modified" impact theory 
comes from the second order expansion of Uj(t). This 
expansion should be valid whenever the quantity 
(i/h) fV(t')dt f is small compared with unity. For com- 
pleted collisions this quantity was estimated (cf. 
eq (1-68)) by (Vr/h) and the expansion was justified 
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except for Strong collisions. For collisions which 
cannot he completed during the time of interest 
(i.e., collisions for which tAco>1), we have (l/h) 
JV(t f )dt — {VjhAo)) and it remains to be seen whether 
this is small or not. This situation may arise in the far 
wings of a line where the electron-atom collisions take 
on a quasi-static aspect (cf. sec. 3.4). In this region, 
the intensity at a point Act) is primarily determined by 
those interactions which produce a Stark splitting 
hAay; that is, by those interactions for which V — hAto 
(see p. 23 of ref. [4]). For these quasi-static interac- 
tions we have (l/h) J V{t')dt' ** (V/hAo)) = 1 and 
the series expansion of Uj(t) breaks down [29]. 
To determine where this breakdown occurs, we use 
t — (po/fav) {vo — v ay for neutral radiators) and rAco > 1 
to obtain p > (t> av /Aco); substituting this result in 
hA(D — V — e 2 /?o/p ( i, we obtain Aco — {hvrJe-Ru) = Aco, 
where Ao> r is the Weisskopf frequency for electrons. 
Using v ax = (SkT/m) 11 ' 2 and ao= (h 2 /me 2 ), we have 
Acor— (3A:77fi) (oo//?o); thus we see that the second 
order expansion of Uj(t) limits the application of the 
impact theory to the region 



< Aw ^ Aav= (3kT/h) (ao/Ro). 



(85) 



In principle this range of applicability could be ex- 
tended by using the full exponential for Uj(t) rather 
than its series expansion. In this regard it should be 
noted that the derivation of the Lewis cutoff is based 
on the second order expansion and, if this expansion 
is to be improved, it may be necessary to improve the 
Lewis cutoff as well. 

7.4. Discussion of the One-Electron Theory 

In the previous section it was shown that the im- 
pact theory is limited to the center region of a line 
profile because of the second order expansion of 
Uj{t). The one-electron theory is designed to be 
valid in the line wings and it is expected that its region 
of validity will overlap with that of the impact theory 
so that the complete line profile may be described 
without ambiguity. To insure that the one-electron 
theory will be valid in the line wings, we used the 
exponential form of Uj(t) which is valid in the wings. 
In the derivation of this exponential form it was 
assumed that both the ion field splitting and the time 
ordering operator could be neglected. 

In section 3.3 it was shown that the time ordering 
operator is not necessary if the second order expansion 
of Uj(t) is justified, and in section 3.4 it was shown 
that time ordering is not necessary for quasi-static 
interactions. It thus appears that the only errors 
resulting from the neglect of time ordering will come 
from the strong collisions. Since the average effect 
of strong collisions is to reduce Uj(t) to zero (see 
sec. 7 of I), the errors made in these collisions should 
have a negligible effect on (Uj{t)). 

To find out when we may neglect the exponential 
exp [iteff(Z„/—Z„/')/h]vihich appears in the operator 
U n (t) used by the one-electron theory, eq (28), we 
consider a cutoff that is similar to the cutoff i;/ow ,,/ 
which approximates this exponential in the impact 



theory (see eq (75)). Using t=l/Aa) and gj„/. „/■(#) 
= e < £{Z„/—Z„//)jh, we note that this exponential is 
essentially unity when (co,,/- ,,/VAaj) < 1 : when 
((*)„/( „/7Ag>) > 1 the exponential oscillates rapidly 
and there is only a negligible contribution from 
U„(t). Since U„(t) appears under a microfield average, 
we will consider a cutoff gf r = ftAa>/e (Z„/-~ Z,j) in 
the range of this microfield average. We replace the 
exponential by unity and average over ion field 
strengths £f<£V; this cutoff replaces the integrand 
by zero when (a»„4 „/'/£&))> 1, thus approximating 
the effect of the rapidly oscillating exponentials, 
and when (a> „/ \ „/, ■jAco) < 1 , the exponentials are 
effectively replaced by unity. This cutoff has no 
effect on the static ion broadening because, in the 
one-electron theory, the 1 line shape is a sum of two 
terms (cf. eq (SI)) and the term which describes the ion 
broadening (the .AftAiolZ „/) term) is not affected by 
the microfield average of the electron broadening 
term (the F\(t) term) where the cutoff is applied. 
If it should happen that^>gf av , we may as well 
ignore the cutoff and average over ^g 7 ^ oc because 
^(g 7 ) will be negligibly small for %> % v . That is. the 
exponential may be replaced by unity for all ^when- 
ever gf r = fiAcoleiZ,,/— Z „/•) >g ? av or, as it was stated 
in eq (29), whenever 



Aoj> eg^\Z„/-Z„t)lh> 



(86) 



It should be noted that, as a practical estimate of 
the region of validity for neglecting the effect of the 
ion field on electron broadening, eq (86), is usually 
too restrictive. For higher lines (large // ) which have 
a large number of Stark components (allowed 
\/*A) —* | 1) transitions) a more realistic estimate 
is obtained by considering the average Stark splitting 



Aco > e%, xs Z,Jfx 



(87) 



which is usually on the order of the halfwidth [30]. 
The use of Z av rather than the maximum value of 
(Z„/— Z,,/-) is more appropriate for the higher lines 
because, for the larger values of Z„/, the quantity 
e% iXX Z„/,lfi corresponds to a very large frequency and 
the electron broadening is negligible for some values 
of Aoj less than eg, XK Z„/ lh (for large Z„/, ^(ftA<o/eZ„/) 
is much larger than the electron broadening term for 
a wide range of Ao> which are less than e% . lXX Z /t /,\h). 
Since we wish to use a classical potential in Uj(t) 
to describe the quasi-static interactions, it is necessary 
to justify the use of the classical wavepackets which 
give rise to a classical potential (in sec. 8.3 of I, we 
justified the use of classical wavepackets only for 
dynamic perturbers). In section 8.3 of I, it was shown 
that we may represent the perturbers by nonover- 
lapping classical wavepackets of width A.r and momen- 
tum uncertainly Ap <^ p if 



^ av < Ax<n 



1/3 



(88) 



Since we have^av^/? l:{ for most problems of inter- 
est, this criterion is easily met. We thus need only 
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show that the radiator and perturber wave functions 
do not overlap during the time of interest. Using the 
notation of section 8.3 in I, this criterion may be written 



[Ajt-h (fA/>//w)+K<>] <p,» 



(89) 



R» 



: //"do — •/ 



l h 2 lme 2 denotes the effective 



extent of the radiator wave function. Using A/; <^ p 
and f <T= (po/vo), we see that (tAp/m) < p {) hence 
the (thp/m) term in eq (89) may be dropped. Using 
V — hAa) (for quasi-static interactions) and V~ e 2 /?o/pii 
we obtain pi =* (e 2 /?o//iAa)) = R%(e 2 l^ 2 ao) (1/ftAcu) 
= /?o (2Ey-J/<-fiA(o) where E yc = e 2 l2ao is the ionization 
energy. Since hko) < E Jn 2 , we see that fj) <^ p<> 
hence the /?<» term in eq (89) may also be dropped 
and we need only show that 



Ay < p„. 



(90) 



Again using the value pi— (e 2 R { JfiAa)) given above, 
as well as A|^= (h 2 l3mkT) and a {) = h 2 lme 2 , we find 
that p?)=A! v (AT/ftAco )(/?o/ <7<>). Since /?<>>«<) and 
ftAoj <^ AT, we see that p () ^>A av hence it is possible to 
find a A.r which will satisfy both eq (90) and eq (88). 
The use of a classical potential function is thus justi- 
fied for quasi-static electron-atom interactions. 
The results of the one-electron theory are based on 
the approximation F(t) — \+ F\(t) which treats the 
electron-atom interactions as though they are inde- 
pendent of one another. In section 4 it was argued 
that this is probably a good approximation, if the 
times of interest are shorter than the mean free time 
between strong collisions. To obtain a validity criterion 
for this approximation one sometimes identifies the 
collision frequency with the halfwidth Aa>i/2 for the 
line under consideration; the requirement that the 
time of interest, 1/Aco, be less than the mean free time 
between collisions is then written [3] 



Aw > Aw 



1/2 



(91) 



Although this argument is somewhat loose, it will be 
shown that eq (91) actually does provide a useful 
statement of the region of validity for the one-electron 
theory. 

To derive a validity criterion for the one-electron 
theroy, we note that the approximate F(t) used in 
this theory, eq (51), is formally the same as an expan- 
sion to lowest order in h\{t) of the approximate 
F(t) used in the impact theory, eq (61). That is, the 
single particle approximation is simply a small hi(t) 
approximation which is valid when n\h\(t)\<\. 
From the definition of h\{t) (see eq (62)) it is clear 
that n\h\{t) \ < 1 for small /. For large t the impact 
theory is valid and we use eqs (66) and (71) to obtain 
n\h\{t) | ^ /0„. Replacing / by the time of interest, 
1/Ao). the validity criterion n\h\(t)\<\ becomes 

Aco ><!>„. (92) 

Thus the one-electron theory should be valid when 
both eqs (86) and (92) are satisfied (or, for higher 
lines, when eqs (87) and (92) are satisfied). 



If the average ion field splitting is small. (e2T av Z*//ft) 
< <!>„, or if the linear Stark effect due to the ion field 
does influence the line shape (e.g.. isolated lines), 
the halfwidth is determined [30] by the matrix ele- 
ments of <P (I (see sec. 4.7 of ref. [2]). In this case, the 
region of validity is given by eq (91) as Aco > Aa>i/2. 

For lines where the average ion field splitting is 
large, (c^ av Z av //? ) > *„ (e.g., higher hydrogen lines) 
the halfwidth is determined primarily by ion broaden- 
ing [30 1 and Aw,/2 - (e%. lXS Z. tW lti) . In this case eq (87) 
is more restrictive than ecj (92) and the region of 
validity is again given by Ao> > Aoh/2. 

The low lying hydrogen lines are determined by a 
small number of Stark components and, for many of 
these lines, a discussion of the halfwidth is of little 
value; lines such as Ly-/3 and H-/3 do not even have a 
well-defined halfwidth. Lines such as Ly-a and H-a 
have a very intense Stark component which is not 
shifted by linear Stark effect and the halfwidths of 
these lines are determined primarily by the electron 
broadening of this unshifted component. In such a 
case, eq (92) gives Ao> > Acoj/2. However, the remain- 
ing Stark components are very important in determin- 
ing the line profile outside of the relatively narrow cen- 
tral component and, for these shifted components, eq 
(86) becomes more restrictive than the condition 
Aw > Aah/2 obtained from eq (92). It should be noted 
that eq (86) is probably not very much more restrictive 
than the condition Aco>Agoi/2; for example, in the 
cases where Ly-a has been studied experimentally 
[31, 32), eq (86) gives approximately Aco>cu,,/10 and 
(Oy,/10 differs very little from Acoi/2. 

From the above discussion we conclude that the 
one-electron theory will be valid in the region 
Ao) > Acl>i/2 for isolated lines and for overlapping lines 
which have a large number of Stark components. 
For overlapping lines with only a few Stark com- 
ponents, the condition Ao> > Acoi/2 may no longer be 
applicable and one should consider eqs (86) and (92) 
for each Stark component. 

8. Comparison of the Impact and 
One-Electron Theories 

The calculations in both the impact and one-electron 
theories are based on binary collisions. The one- 
electron theory treats the electron-atom interactions 
as though their effects are additive and for this reason 
the one-electron theory is limited to the line wings 
(short times of interest) whereas the impact theory 
is not. In principle, the impact theory (with a Lewis 
cutoff) should be valid over most of the line profile; 
however, in practice it is usually necessary to expand 
the time development operator, Uj(t) , in order to 
make practical calculations, and it is this expansion 
which limits the impact theory to the center region of 
a line profile. 

It was shown that the series expansion of Uj(t) 
could not adequately describe either the strong col- 
lisions, (Vr/h) > 1, or the quasi-static interactions, 
(i/h)fV(t)dt — 1, which are important in the line 
wings. In our version of the one-electron theory this 
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difficulty is avoided through the use of an exponential 
expression for Uj(t) (previous versions of the one- 
electron theory have used either a completed collision 
assumption [5] or a series expansion [3] of Uj(t) 
hence they cannot be valid in the quasi-static wings). 
While this exponential form is not rigorously correct 
for strong collisions it does have the following useful 
properties: 

(1) It is correct for weak collisions, (Vr/h) < 1, 
since it reproduces the series expansion for such inter- 
actions. 

(2) It is correct for quasi-static interactions whereas 
the series expansion is not. 

(3) For strong collisions, (Vr/h) > 1, the exponen- 
tial would give a much more realistic approximation 
to the exact Uj(t) than the series expansion because 
the exponential retains unitarity (i.e., no transition 
probability can exceed unity). Thus, if one chose to 
use classical path methods for calculating Uj{t) , the 
strong collision cutoff discussed in section 7 of I, 
would be essentially reproduced (i.e., the exponential 
would oscillate rapidly and contribute essentially 
nothing for strong collisions). The errors resulting 
from this treatment of strong collisions should be 
negligible because the averaged Uj(t) for strong col- 
lisions is essentially zero. 

It is interesting to note that if we use a second order 
expansion of Uj(t) in the one-electron theory, we ob- 
tain the asymptotic wing expansion which is used in 
deriving the Lewis cutoff [21]. That is, the Lewis 
cutoff, which is applied to the impact theory in an ad 
hoc manner, is an integral part of the one-electron 
theory. 

Since the Weisskopf frequency, Aco c , for electrons 
is generally much greater than the halfwidth, there 
should be a broad range of overlap between the re- 
gions of validity 0^Aco<Aco r and Acoi/ 2 <Ao) for the 
impact and one-electron theories. Nevertheless, in 
our derivation of the impact and one-electron theories 
we have tried to stress the formal connection between 
the two in order to provide a better understanding of 
their similarities as well as their differences. We feel 
that this is particularly important since it would be 
desirable to formulate a classical path theory which is 
valid in the line center as the impact theory is, and 
which makes a smooth transition to the line wings as 
the one-electron theory does. It is hoped that the 
present approach will make it possible to remove some 
of the approximations which currently limit these 
classical path theories. 

9. Appendix. Collision Variables 

We wish to use the natural collision coordinates in 
evaluating the integral 

Fi(t) =n \dy x W{vi) J dxj 

[expf-(i/ft) rPFi(R,x, + v,OA']-ll (A. 
where 



K,(R,x, + v,*') = e 2 R-(xi + v,t)/|xi + v,t| 



(A.2) 



1) 



For fixed values of xi and vi, we note that the head of 
the vector (xi+vi*') traces out a straight line if/' is 
varied from — oo to + °°. We shall call this lineL(xi, vi). 
The integral over t' , in eq (A.l), will cover a segment of 
this line. For fixed xi and Vi, we can define a vector 
po from the origin of coordinates to the closest point 
on L; this vector will be perpendicular to vi. We can 
also define a scalar to by 



xi = po + Vifo. 



(A.3) 



For a fixed vi, we see that, as xi varies over all space, 
|vi|fo varies on the range (— °°, + 00 ) and po=\po\ 
varies on the range (0, oo). We may thus change vari- 
ables from xi to the cylindrical coordinates (po, £o, «) 
where vi defines the polar direction and a is an angle 
which ranges from to 27T. The Jacobian for this trans- 
formation is such that dxi — » | V\ \ podpodtoda. In these 
coordinates the potential function becomes 



^i(R,po,Vi,*' + to) = e 2 R 

•[po + Viit' + hWiri + W + to) 1 ]* 1 *. 



(A.4) 



We may express d\\ in spherical coordinates, vldvo 
sin /3 df3 dy, where the angles f3 and y denote the ori- 
entation of vi relative to the direction of R. Since po 
is defined to be orthogonal to vi, the integrals over the 
Euler angles a, /3, and y average over all possible 
orientations of the (po, Vi) frame of reference (collision 
axes) relative to the direction of R. Equation (A.l) may 
now be written 

Fi{t)=n dCl vW(vo)dvo podpo du 

J JO J J-oo 

[exp{-(i/ft)J Wi(R, p* f v ,*' + *)*'} -l| (A.5) 

where dil = sin (3 dad/3 dy. 

It is clear that the line L represents the trajectory 
of a perturbing electron, vo is the speed of the electron, 
po is the impact parameter, and Jo is some starting time 
for the collision (or, if the t' limits of integration are 
(—oo, -foo) in a completed collision, t may be regarded 
as the time of closest approach). The orthogonal 
"collision axes" to which this collision is referenced 
are po and Vi. The Euler angles Ci denote the orienta- 
tion of the collision plane relative to the direction of 
R; the integral over ft may therefore be interpreted as 
an average over all possible orientations of the electron 
trajectory. 
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